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Abstract. Adaptive optics observations of NGC 1068 
have allowed to reach with a 4 m class telescope, 
diffraction-limited images at 3.5 /xm (FWHM=0.24") and 
4.8 /Ltm (FWHM=0.33") of the central arcsec region. 
These observations reveal the presence of an unresolved 
core (radius less than 8 pc at Half Maximum) and an 80 pc 
sized "disc-like" structure at P.A. '^100°, which is inter- 
preted as the dusty/molecular torus invoked in the AGN 
unification scheme. They show as well an extended emis- 
sion region along the NNE-SSW direction, of 100 pc full 
size, most probably associated with dust in the NLR. The 
position of the unresolved core at 3.5 and 4.8 /im is found 
to be coincident with that of the core observed at 2.2 
/xm and outlines the location of the central engine in the 
AGN of NGC 1068. Considering as well previous AG ob- 
servations at 2.2 /Ltm we infer that there must exist a very 
steep gradient of the dust grain temperature, close to the 
central engine. At a distance of 30 pc from the central 
heating source, the dust grain temperature deduced from 
the [L-M] color is found to be ~500K. The mass of warm 
dust {Tgr >500K) in the 0.6" diameter core is found to 
be ~0.5 Mq. The spectral energy distribution from 1 to 
10 /im is provided for the 0.6" diameter core. These re- 
sults are briefly discussed in the context of current torus 
models. 
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1. Introduction 

Series of observational facts assembled over the past 
decade on Active Galactic Nuclei have led to the so- 
called "unified" model of AGN (for a review of these facts, 
as well as for the detailed characteristics of the unified 
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model see Krolik 1999). Our current specific interest in 
the unified model is that the central engine (black hole 
and accretion disk) and its close environment (dense gas 
clouds emitting the broad lines which constitute the broad 
line region, BLR) are embedded within an optically thick 
dusty/molecular torus. Along some lines of sight, the torus 
obscures and even fully hides the central engine and the 
BLR. 

In that respect, the case of NGC 1068, a bright Seyfert 
2 active galaxy, is particularly enlightening. The spec- 
trophotometry in polarized light reveals the presence of a 
hidden BLR (Antonucci & Miller 1985), the conical shape 
of the narrow line region (NLR) ~ both on large (Pogge 
1988) and small (Evans et al. 1991) scales - indicates that 
the ionizing radiation is coUimated by an opaque blocking 
torus and, finally, the symmetry center of the UV/ optical 
polarization map (Capetti et al. 1995) is found to be co- 
incident with the radio core (Gallimore et al. 1996a), the 
12.4 fim peak (Braatz et al. 1993) and the maser emis- 
sion (Gallimore et al. 1996b), suggesting that this is the 
location of the hidden true nucleus. 

This object appears then particularly suitable for un- 
veiling the putative torus through its infrared emission, 
which, according to current models should be quite strong 
{e.g. Pier & Krolik 1992a, b, 1993; Granato & Danese 
1994; Efstathiou, Hough & Young 1995; Granato, Danese 
& Franceschini 1997). In this search, high spatial resolu- 
tion is a requisite in order to locate very precisely and 
to characterize the structure of emission sources. Hence, 
adaptive optics (hereafter abbreviated AO) in the 1-5 
/xm window is the tool. At the distance of NGC 1068 
(14.4 Mpc), 1" is equivalent to 72 pc (assuming Ho=75 
km/s/Mpc), allowing to reach a spatial resolution of a 
few parsecs. Through AO observations, simultaneously in 
the visible and near-infrared, the 2.2 ^m peak emission has 
already been found to be offset by ~0.3" S of the optical 
continuum peak as defined by Lynds et al. (1991) and to 
be coincident with the previously identified "hidden true 
nucleus" (Marco, Alloin & Beuzit 1997). In the current 
study, we are presenting new results obtained at 3.5 and 
4.8 fim with ADONIS, the AO system working at the ESO 
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Table 1. Summary of the data sets 



filter 


A 


Tint 


Ttotal 


Strehl 


visible 


airmass 




[Mm] 
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seeing ["] 
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3.48 


6.0 
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28 
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1.23 
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4.83 


1.5 


720 
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0.55 


1.15 


PAH 


3.31 


60 


1440 


25 


0.66 


1.25 



3.6 m telescope on La Silla and fully described in Beuzit 
et al. (1994). 

2. Observations and data reduction 

The observing run took place from august 13 to 19, 1996, 
under excellent seeing and transparency conditions (Table 

0). 

The AO correction was performed on the brightest 
spot of NGC 1068 in the visible continuum light (Lynds 
et al. 1991). The wavefront sensor (EBCCD after a red 
dichroic splitter) has a pixel size of 0.7" and takes into ac- 
count the gravity center of the light within a 6" diameter 
circular entrance. Due to the pixel size of the wavefront 
sensor, the contribution of the continuum and the lines 
from the central 50 pc around the Lynds et al. peak (1991) 
both fall in one pixel: so the contrast is maximum. 

The detector was the COMIC camera (Lacombe et al. 
1997), at the f/45 Cassegrain focus, which provides an 
image scale of 0.1"/pixel, resulting in a field of view of 
12.8"x 12.8". 

NGC 1068 was observed in an imaging mode through 
the standard spectral L (3.48 nm), L' (3.81 /im), and M 
(4.83 fim) bands and through a circular variable filter for 
the PAH (line 3.3 /im rest wavelength, and continuum). 
Through the L, L' and M bands, we observed in a chopping 
mode, alternating object and sky images by the use of a 
field selection mirror. We chose an offset of 10" to the N 
and 10" to the W. 

During the six nights, the visible seeing was measured 
by the ESO differential image motion monitor. It was 
excellent during four nights, ranging from 0.4" to 0.7". 
Therefore, the efficiency of the AO correction was optimal 
and the images obtained with COMIC were diffraction- 
limited. However, the gain of the intensifier was not opti- 
mized at that time, which resulted in Strehl ratios lower 
than expected. 

In order to minimize position offsets between the cal- 
ibration star and the AGN, we selected a reference star 
within 2 degrees of the target. For both the galaxy and 
the point spread function (PSF) reference star, the air- 
mass was at most of 1.3, ensuring differential refraction 
effects to be negligible (less than one pixel). 

Individual exposure times were chosen so as to observe 
under conditions of background limiting performances 
(BLIP). In this way, the readout-noise is dominated by 



the background photon noise, and we just take an aver- 
age of the images to improve the signal-to-noise ratio. We 
observed several photometric standard stars to obtain a 
precise flux calibration and another star to determine the 
PSF for later deconvolution. 

Standard infrared data reduction procedures were ap- 
plied to each individual frame, for both the galaxy and the 
reference stars : dead pixel removal, sky subtraction, flat 
fielding from sky images at each wavelength. As the AO 
system is compensates for image shifting, no additional 
shift correction was applied. 

Thanks to the length of our observing run, we obtained 
a largely redundant data set. We discovered in particular 
that we had experienced a problem of astigmatism during 
2 nights, due to an unappropriate tuning of the visible 
wavefront sensor leading to a very low signal to noise ratio 
on the Shack-Hartmann analysor (Alloin & Marco 1997). 
Being aware of this flaw, we decided to review in depth the 
AO optimization parameters for the entire data set and 
to remove all suspicious blocks. In addition, we applied a 
selection procedure of 32-images data cubes, based upon 
the seeing value and the Strehl ratio. The corresponding 
equivalent integration times are reported in Table |l|. 

3. Wavefront analysis sensing 

To better interpret the NGC 1068 observations with ADO- 
NIS, we need to track precisely which visible image of 
NGC 1068 is seen by the wavefront analysis sensor (WFAS) 
in its evaluation of the AO correction. 

In order to recover this information, we have used HST 
images with 45 mas pixel size both in the continuum and 
in the lines (F791W, F547M F658N, and F502N), prop- 
erly aligned and flux scaled to a same 1 sec exposure at 
the entrance of the telescope, to construct a composite 
image which mimics the one seen by the WFAS. The im- 
age alignment was performed using point-like sources in 
the field of view, which allows a registration within 10 mas 
(Tsetanov, private communication). The flux scaling has 
included corrections for the exposure time of each image 
and the mean HST camera efficiency over the correspond- 
ing filter. Then, each of these corrected images has been 
weighted by the corresponding mean WFAS wavelength 
response, before addition. We find indeed that the light 
received from NGC 1068 on the WFAS is largely domi- 
nated by continuum light and is not much weighted by 
the [O HI] - line light distribution. 

And finally, to mimic the image seen on the WFAS 
CCD, from which the AO corrections are computed, we 
have degraded the HST composite image by a "seeing" 
effect of '--^1": this includes both an atmospheric effect 
(~0.6" according to the mean seeing value -FWHM- mea- 
sured on the selected nights of the observing run) and 
the instrumental spread function (0.85"), which is rather 
large because of the photocathode -CCD spacing. As the 
WFAS takes into account, to compute the AO correction, 
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Fig. 1. Left: L band image of NGC 1068. Right: M band image of NGC 1068. North is up, East is to the left. 



the light centroid over a 6" circular diaphragm, we have 
derived the hght centroid within a 6" aperture on the de- 
graded composite HST image. This position is the refer- 
ence for the WFAS. 

We need also to locate the centroid of the degraded 
HST composite image with respect to the Lynds visible 
peak, as measured from the HST image F547M. We find 
a slight offset between the visible centroid and the Lynds 
visible peak: the visible centroid (using NGC 1068 for the 
correction) is located 99 mas to the N and 86 mas to the E 
of the Lynds peak. This offset is taken into account in our 
estimate of the NGC 1068 infrared sources location (Sect. 
4.1). The offset may differ from one AO system to another. 
In the case of ADONIS, the offset is largely due to the 6" 
entrance to the WFAS. 

4. Data analysis 

We present in Fig. |l| the L and M band images of 
NGC 1068. We used a magnitude (log) scale because of the 
high dynamics of the images provided by the AO. The im- 
ages have been deconvolved using a Lucy-Richardson algo- 
rithm (MIDAS package). We have also observed ngc 1068 
in the L' band: this image is very similar to the one ob- 
tained in the L band. Thus, our current data analysis and 
subsequent discussion will be based on the two L and M 
bands, only. 

The L and M band images show: 

i) an unresolved core down to the resolution (FWHM) 
of 0.24" (16 pc) and 0.33" (22 pc), respectively at 3.5 and 
4.8 /im. This core has already been observed at 3.6 /im by 
Chelh et al. (1987) and at 2.2 fj.m by Marco et al. (1997), 



Thatte et al. (1997) and Rouan et al. (1998). The latter 
give an upper limit of the core size (FWHM) of 0.12" (less 
than 8 pc). 

ii) an elongated structure at P.A. ^100° particularly 
prominent in the M band, but also quite well outlined in 
the L band. This structure is obviously coincident with the 
structure seen in the K band by Rouan et al. (1998) and 
is roughly perpendicular to the axis of the inner ionizing 
cone (P. A. =15°, Evans et al. 1991). It extends in total 
over ~80pc, with a bright spot at each of the ^E and 
^W edges at a radius of ~25pc from the central engine. 

iii) an extended emission along the NS direction, al- 
most aligned with the radio axis and the ionizing cone 
axis. At low level isophotes (in particular in the L band), 
a change in the direction of the axis of this emission can 
be noticed, reminiscent of a similar change of direction of 
the radio jet (Gallimore et al. 1996a). 

Down to faint levels, the 4.8 /im thermal infrared 
emission appears to be extended over ^3" in diameter 
(-210pc). 

It is striking that the two different AO systems used, 
ADONIS for the L and M bands and PUEO for the K 
band, reveal a similar structure of the AGN environ- 
ment. These AO systems are using WFS of different 
types, Shack-Hartmann for ADONIS and curvature for 
PUEO, as well as deformable mirrors of different types, 
piezo-stack for ADONIS and bimorph for PUEO. The 
Lucy-Richardson deconvolution applied on both data 
sets uses PSFs obtained in two different ways: we used 
an observed stellar PSF in the case of the ADONIS 
data set - as the L and M band data are less sensitive 
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Table 2. Photometric data for NGC 1068. Observed fluxes 
are in Jy. 
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aperture radius 






0.3" 


0.4" 0.5" 0.6" 0.9" 


1.5" 


L 


0.92 


1.34 1.78 2.09 2.76 


3.40 


M 


2.27 


3.23 4.41 5.30 7.18 


9.73 



to rapid PSF fluctuations - and we used the PSF 
recovered from the AO loop parameters in the case of 
the PUEO data set. The two AO experiments differ 
in many aspects, while leading to a similar result for 
the structure of the AGN dusty environment. Therefore 
we are quite confident that this structure is real and 
not hampered by significant AO artifacts (Chapman 
et al. 1999). Finally it should be noticed that the high 
resolution image of the AGN in NGC 1068 obtained in 
the K band with the AO system at the Keck telescope 
(www2. keck, hawaii. edu/realpublic/ao/ngcl 068. html) re- 
veals a comparable structure, pending that a precise 
orientation and a scale be provided for the Keck AO data 
set. 

^.1. Location of the emission peaks at 3.5 and ^.8 iiva 
and nature of the unresolved core 

As it was not possible to observe simultaneously in the 
visible and in the infrared, we took advantage of a char- 
acteristic feature of AO systems which is to preserve the 
optical center for all objects: the infrared camera field has 
a position fixed in regard to the centroid of the visible 
counterpart of the object observed. Indeed, by observing 
a star (PSF or photometric standard), we determine a 
reference position in the infrared image to within the pre- 
cision we are aiming at in this study (better than one 
infrared camera pixel, 0.1"). Any offset of the galaxy in- 
frared peak relatively to the star infrared peak would then 
reveal an intrinsic offset between the galaxy infrared light 
peak and the galaxy visible light peak. This is a method 
for positioning infrared versus visible sources in the AGN. 

To improve the precision, we have fitted the PSF and 
the NGC 1068 emission peaks by Gaussian profiles. The L 
and M band peaks in NGC 1068 are coincident within the 
positional precision given above. We have also derived the 
position of this L and M peak in NGC 1068 with respect 
to the visible peak, following the procedure described in 
Sect. 3: it is offset by 0.3 ± 0.05" S and 0.1 ± 0.05" W of 
the visible continuum peak and therefore is found to be 
coincident with the K band emission peak (Marco et al. 
1997), within the error bars. 

Therefore, the compact core at 3.5 and 4.8 /^m can 
be identified with the unresolved core detected at 2.2 ^m 
(Marco et al. 1997; Thatte et al. 1997; Rouan et al. 1998), 
itself found to be coincident with the mid-infrared emis- 
sion peak at 12.4 /xm (Braatz et al. 1993), the radio source 



SI (Gallimore et al. 1996a, b) and the center of symmetry 
of the UV polarization map (Capetti et al. 1995). This 
strengthens considerably the interpretation of the core in- 
frared emission originating from hot/warm dust in the im- 
mediate surroundings of the central engine. 

^.2. The torus-like structure at 3.5 and 4-8 /ini 

The location, position angle and extension of the 
P.A. ^100° structure are strongly suggestive of a 
dusty/molecular torus. The two bright spots on the edges 
of the structure outline the "disky" nature of the torus, 
up to a radius of ~40pc from the central engine. This 
dusty/molecular torus would be responsible for the coUi- 
mation of UV radiation from the central engine, leading 
to the ionizing cone (Pogge 1988; Evans et al. 1991). The 
overall spatial extension of the torus is found to be ~80 pc 
at 3.5 and 4.8 /im, while it appears to be slightly smaller at 
2.2 fiiri, ^50 pc. Under the very simple assumption of opti- 
cally thick grey-body dust radiation (Barvainis 1987), it is 
well understood that the emission at 2.2 fim traces hotter 
dust (T~1300K) than the emission at 4.8 fim (T-600K). 
The observed difference in size would then signal the ex- 
istence of a temperature gradient of the grains across the 
torus. 

4.3. The North South extended emission 

The NS extended emission (overall extent ^^3" down to 
faint emission levels) is also detected at 2.2 /im on a similar 
scale (Rouan et al. 1998) and at 10 and 20 /im on a larger 
scale, although along a similar P.A. (AUoin et al. 1999) . 
This structure can be related to the emission of hot /warm 
dust associated with NLR clouds identified in the north- 
ern side of the ionization cone from HST data (Evans et 
al. 1991) and hidden behind the disc of the galaxy in its 
southern side. Additional local heating processes, e.g. re- 
lated to shocks induced by the jet propagation , might be 
at work as well along the NS extension. The latter sugges- 
tion stems from the conspicuous change of direction of the 
emission at 3.5 /im, following that of the radio jet. Indeed, 
Kriss ct al. (1992) have shown through the analysis of line 
emission that in NGC 1068 the emitting gas in the NLR is 
partly excited through shocks triggered by the radio jet. 

5. Fluxes, SED and variability 

The spectral energy distribution (SED) of the central re- 
gion of the AGN is an essential parameter in the model- 
ing. To derive this quantity, spatial resolution is obviously 
needed to disentangle the different sources of emission - 
dust, stars, non-thermal source.. - and, in that respect, 
AO observations bring precious informations. 

Fluxes at 3.5 and 4.8 /xm have been measured through 
circular apertures centered on the near-infrared peak, with 
a radius varying from 0.3" to 1.5" (22 to 100 pc). They are 
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depicted in Table g. A weak PAH line emission at 3.3 /im 
has been detected as well, although no flux calibration is 
available for this observation, unfortunately . 

The aperture flux density as a function of radius, over 
the region 22 pc < r < 100 pc, can be fitted with a power 
law: we find Fl oc r^^"^ and Fm oc r"^°*' (where the flux 
unit is Jy/arcsec^). 
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Fig. 2. SED of the central 0.6" diameter core (the value 
for the N band is given for 0.8", hence the dashed line 
between the M and N bands). The dotted line represents 
an interpolation for establishing the lower-limit of the dust 
contribution in the K band. 

From the set of AO images (this paper and Rouan 
et al. 1998), as well as high resolution images obtained at 
10 /j,m (Alloin et al. 1999), we have reconstructed the SED 
of the core emission through an 0.6" diameter diaphragm, 
as shown in Fig. g. Therefore, this SED refers only to the 
central engine and its immediate environment (including 
the inner parts of the dusty/molecular torus as well as 
some contribution from the NS extended structure). In 
this plot, the contribution from the stellar component has 
been removed in the J and H bands, following the spatial 
profile analysis performed by Rouan et al. (1998), while 
this is not the case in the K band. The part of the near- 
infrared to mid-infrared emission which is arising from 
hot to warm dust grains can be represented by a series 
of grey-bodies of different temperatures and this summa- 
tion is expected to result in a smooth distribution. Yet, 
an emission bump at 2.2 /im is observed, which could be 
interpreted as the unremoved stellar contribution in the 
K band, within the 0.6" diameter diaphragm. Under this 
assumption, an upper limit of 75% for the stellar contri- 
bution to the flux at 2.2 /xm, can be derived. This upper 
limit remains far larger than the 6% effectively derived by 



Thatte et al. (1997) within a 1" diameter diaphragm, on 
the basis of the dilution of the equivalent width of a CO 
absorption feature arising from cold stars. 

In any case, a revision of the AGN modeling for 
NGC 1068 should incorporate the 1 to 10 ^m SED derived 
for the innermost region around the central engine and 
therefore less affected by dilution from other surrounding 
components (in particular cleaned from the stellar contri- 
bution) . 

It is interesting as well to compare these 1996/1997 
measurements in the infrared to those performed by Rieke 
& Low as early as 1975. Therefore, we have derived, 
for their 3" diameter diaphragm, the 1996/1997 observed 
fluxes at 1.25, 1.65 (both uncorrected for the stellar contri- 
bution), 2.2, 3.4, 4.8 and 10 yum, from Rouan et al. (1998), 
from the current data set and from Alloin et al. (1999). 

Before examining the temporal behavior of the near- 
infrared emission in NGC 1068, the 1997 flux measure- 
ments in the K band can be compared to previous determi- 
nations. We consider again a 3" diameter diaphragm and 
compare the 1997 AO measurement to classical photomet- 
ric measurements. Generally those are performed through 
much larger diaphragms. However Penston et al. (1974) 
provide a full set of measurements through diaphragms 
ranging from 12" to 2" diameter. Leaving aside the mea- 
surements from this paper which have been flagged down 
for bad transparency or other flaws, we obtain the K mag- 
nitude offset when one moves from a 12" diaphragm to a 
3" diaphragm, AK = 0.85. This magnitude offset is not 
expected to vary with time because it refers to the outer 
parts of the AGN (r > 110 pc). Then, from the 12" di- 
ameter diaphragm measurements by Glass (1995), who 
analyzed the variability properties of NGC 1068, we can 
infer/extrapolate the K magnitude at the date of the AO 
measurement by Rouan et al. (1998): K = 6.81. Applying 
the magnitude offset computed above between the 12" and 
3" diameter diaphragms, we predict that the K magnitude 
should be of 7.66 at the date of the Rouan et al. (1998) 
observation, while the K magnitude measured is of 7.26. 
This agreement is quite satisfactory, given the rather large 
error-bars involved in the Penston et al. (1974) data set 
which was obtained more than 25 years ago. 

A comparison of the fluxes within a 3" diameter di- 
aphragm at both epochs, 1975 (Rieke & Low 1975) and 
1997 (this paper) is depicted in Table 0. One notices im- 
mediately that a flux increase has occurred over this time 
interval: by a factor 2 between 2.2 to 4.8 /zm, and by a 
factor 1.2, at 10 /im. According to the independent pho- 
tometric monitoring by Glass (1997), the L band (3.5 /im) 
flux has doubled in 18 years, from 1974 to 1992: our result 
is in very good agreement with his finding. 

Because of the possibly complex way through which 
the UV-optical photons illuminate and heat the dust 
grains (direct and/or indirect illumination via scattering 
on the mirror to the N of the central engine or on dusty 
regions further out. Miller et al. 1991), it is not possible to 



O. Marco & D. AUoin: Dusty/Molecular Torus in NGC 1068 



Table 3. Comparison of the infrared fluxes (given in Jy) 
at two epochs in a 3" diameter diaphragm 
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infer the size of the dust component from hght-echo effects 
acting differentiaUy in the near-infrared and mid-infrared 
bands. 



6. Colors, dust temperature and mass of dust 

Following Barvainis (1987), we assume the emissivity of 
the grains to depend on wavelength as A~^'^ and we use a 
simple model of optically thick grey-body dust emission: 

In the general situation of an AGN, the grain tem- 
perature varies strongly with distance to the central en- 
gine, following a power law (Barvainis 1987): Tg^{r) — 



1650L' 



0.18 
UVA6' 



-0.36 



K, where Ljjy^Q is the UV luminosity 



in units of 10^^ ergs s ^ and r the radial distance in par- 
sees. 

6.1. Color gradients in the near-infrared 

6.1.1. Colors of the core 

Owing to our limitation in spatial resolution in the M 
band (FWHM ~ 0.33"), we can measure at best the [L- 
M] and [K-L] colors of the core through an 0.6" diameter 
diaphragm centered on the near-infrared emission peak. 
We find [L-M] = 1.6 ± 0.4 and [K-L] = 1.8 ± 0.2. It must 
be noticed however that the contribution of stellar light 
in the K band (Thatte et al. 1997) has not been removed 
at this stage and that the observed [K-L] color does not 
relate only to the dust emission. 

6.1.2. Colors of the extended structure 

As for the colors of the sources forming the extended struc- 
tures, and again because of different spatial resolutions in 
the K, L and M bands, we have considered the mean col- 
ors in a ring which extends from r = 0.3" to r = 0.5". 
This ring does include the emitting regions forming the 
extended structures along the two directions P.A. ^^100° 
and NS, at 3.5 and 4.8 /im, but excludes in part the sec- 
ondary peaks which delineate the extended structures at 
2.2 ^m. Therefore, at 2.2 /im, the ring includes more of 
the diffuse contribution possibly related with the stellar 
core analyzed by Thatte et al. (1997). The colors found 
for the ring, representative of a mean 0.4" radius, are [L- 
M] = 1.6 ± 0.4 and [K-L] = 2.8 ± 0.2. 



6.2. Dust temperature 

The advantage of interpreting the [L-M] color is that the L 
and M band flux contributions are known to arise almost 
entirely from the dust component. Within the limitation 
in spatial resolution from the M band data set, we do 
not detect any [L-M] color gradient within the central 1" 
region of NGC 1068. 

Under the simple assumption of optically thick grey- 
body emission from the dust component, the observed [L- 
M] color corresponds to a grain temperature Tgr^^SOK. 
The foreground extinction to the core has been calculated 
by several authors (Bailey et al. 1988; Bridger et al. 1994; 
Efstathiou et al. 1995; Young et al. 1995; Veilleux et al. 
1997; Glass 1997; Thatte et al. 1997; Rouan et al. 1998) 
and an estimate of Ay^SO mag is retained. Applying a 
correction for such an extinction, we deduce an intrinsic 
color [L-M] = 0.8 ± 0.4 and Tg^-700K. The absolute lu- 
minosity Ljjv.opt of the AGN in NGC 1068 - assumed here 
to be the unique heating source of the dust grains in the 
AGN environment - can be approached only indirectly 
and is still quite uncertain. From the analysis by Pier et 
al. (1994) who have examined various methods for deriv- 
ing the absolute luminosity of the AGN in NGC 1068 and 
summarize the current knowledge on this question, we de- 
duce Ljjv.opt = 4 10'*'* erg s^*. However it should be noted 
that this figure has been obtained assuming a reflected 
light fraction irefi ^0.01, while the consideration of ion- 
ized gas in regions further out inside the ionization cone 
(ENLR) leads to a value irefi ^0.001 (Bland-Hawthorn et 
al. 1991). Then a value as high as LiJv,opt = 4 10^^ erg s~* 
should be envisaged as well. In addition, most of these es- 
timates have been derived without taking into account the 
fraction of UV-optical flux which provides the dust grains 
heating: already the K magnitude of the innermost core 
(FWHM — 0.12"), 9.3, corresponds to an energy output 
of ^8 10^* erg s^*. It might be important to consider the 
energy radiated in the near- to mid-infrared bands for the 
evaluation of l^uv,opt- In conclusion, the figures currently 
available for Luv,opt in NGC 1068 might be lower limits. 

Still, under the simple model of optically thick grey- 
body dust emission, reaching Tgr^4:80K at r = 28 pc 
requires Ljjv.opt — 8 10^^ erg s~*, a value roughly con- 
sistent with the highest figure given above for l^ijv.opt in 
NGC 1068. This figure goes up to Luv,opt = 6.5 10'*'' erg 
s~^ if the grain temperature is of 700 K at r = 28 pc 
(extinction-corrected estimate), pushing NGC 1068 to the 
limit between AGN and quasars. This question certainly 
deserves further attention and above all the consideration 
of a more elaborated model of the dust region, with regard 
to its geometry and heating. This is beyond the scope of 
the current paper and will be discussed in the future. 

The [K-L] color in the extended structures can be 
contaminated by some stellar contribution in the 2.2 /im 
band. From grains at Tgr^^SOK., which are the dominant 
contributors, we expect a [K-L]rft,si color of 4.0. Given the 
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observed [K-L] value, we deduce that the pourcentage of 
the flux at 2.2 /Ltm which arises from the dust component 
(with at most rgr~480 K) is of 30%. For this estimate, we 
have not considered any correction for extinction. 

How can the lack of [L-M] color gradient between the 
0.6" diameter core and the extended structures be ex- 
plained? Given the unresolved and intense core emission 
at 2.2 ^m (FWHM = 0.12" from Rouan et al. 1998) it can 
be inferred that the hottest dust grains are extremely con- 
fined and located at a radius less than 4 pc. With Luv,opt 
= 8 10'*^ erg s~^, they would be present only up to r = 
1.1 pc. Hence, there must exist a very steep dust temper- 
ature gradient close to the central heating source. Such a 
few parsec scale corresponds to a resolution which is well 
beyond that accessible at 3.5 & 4.8 /im. In fact the L and 
M emission we are measuring in an 0.6" diameter core is 
already strongly dominated by the warm grains. Because 
of this suspected steep temperature gradient, the proce- 
dure applied previously to derive the stellar contribution 
in the 0.3" to 0.5" radius ring cannot be used in the core. 

6. 3. Mass of the hot dust 

The mass of dust associated with the near-infrared emis- 
sion can be estimated only in a rough way, as it de- 
pends on the (unknown) grain composition and grain 
size distribution. Assuming graphite grains and follow- 
ing Barvainis (1987), the infrared spectral luminosity 
of an individual graphite grain is given by: i^';^ = 
4:Tra'^TTQ,jB,j{Tg^) ergs s~^ Hz~^ where a is the grain ra- 
dius, Q,y ~ q„ v'^ is the absorption efflciency of the grains, 
and B,y{Tg,:) is the Planck function for a grain tempera- 
ture Tgr. Following Barvainis (1987), we take a=0.05 /im 
and in the near-infrared, qii — 1.4 x lO^^'', 7 = 1.6 leading 
to Q^=0.058 (for the K band). 

Because no [L-M] color gradient is detected towards 
the 0.6" diameter core, we consider the simple case of 
2 populations of dust grains in that region, hot grains 
at T = 1500K and warm grains at T = 500K, matching 
the extreme values in that region. Solving the equation 
-Fiv.measured = a;F^,i5ooK + yF;,,5ooK for the three bands 
available, K, L and M, one derives 2 10"*^ and 9 10^^, for 
the number of grains at 1500 K and 500 K respectively, in 
the 0.6" diameter core. This indicates that there are '^450 
times more warm dust grains than hot dust grains in the 
0.6" diameter core. The warm dust grains dominate the [L- 
M] color. With a grain density p=2.26 gcm~^, the mass of 
warm dust grains is found to be Af (warm dust)~0.5 Mq. 
This mass is above that of hot dust grains detected in two 
Seyfert 1 nuclei: 0.05 Mq in the case of NGC 7469 (Marco 
et al. 1998) and 0.02 Mq in the case of FairaU 9 (Clavel et 
al. 1989). This result supports the fact that only a small 
fraction of the dust present in the torus is heated close to 
its sublimation temperature. 



7. Grain number density in the warm dust 

In NGC 1068, like in other AGN (see for instance 
NGC 7469), dust is present also in the NLR region (see Fig. 
|l|). It has been shown in Sect. 5 that the infrared emission 
in the L & M bands over the region 22 < r < 100 pc 
follows a power-law. In addition the temperature of the 
grains deduced from the [L-M] color has been found to 
remain roughly constant with radius. Thus, we can derive 
the radial profile of the warm grain number density in the 
NLR region, ri{r) oc r~^ . Under optically thin conditions, 
probably applicable in the NLR, the brightness is directly 
proportional to the grain number density. Then the ob- 
served brightness power-law leads to (3—1.0, suggesting a 
concentrated grain distribution (a uniform grain density 
would correspond to (3=0). In NGC 7469, the warm dust 
component led to a figure for (3 around 1.5 (Marco et al. 
1998). 

8. Comparison with thick tori model predictions 

Several torus models have been developed so far to ex- 
plain the obscuration of the BLR and UV/X-ray contin- 
uum sources along some lines of sight (AGN unification 
scheme). Some of these models are generic, while others 
have been designed to match the case of NGC 1068. 

Pier & Krolik (1992a, b, 1993) propose a thick, parsec- 
scale, uniform density torus illuminated by a nuclear 
source. The dust can be heated up to the effective tem- 
perature of the nuclear radiation at the inner edge of the 
torus. They investigate models with effective temperatures 
between 500 K and 2000 K. Such a model can explain the 
unresolved core observed at 2.2, 3.5 and 4.8 /Ltm with AO 
in the particular case of NGC 1068. Does it explain the 
extended near infrared emission also revealed by these ob- 
servations? Indeed, extended emission over 1" to 2" could 
result from reflected radiation from the torus and/or dust 
in the NLR. Therefore, this model accounts for most of 
the features unveiled with high resolution imaging in the 
near infrared. 

Efstathiou & Rowan-Robinson (1995) propose a model 
with a very thick tapered disk. They assume the melting 
temperature of all dust grains to be identical (1000 K), 
but consider a radial distribution of the grain physical pa- 
rameters (size and chemical composition). In the case of 
NGC 1068, Efstathiou et al. (1995) have shown that the 
torus emission alone cannot account for the total infrared 
emission. They attribute the excess infrared emission to a 
distribution of optically thin dust with /3=2 in the NLR 
region. Their model is in disagreement with the steep grain 
temperature gradient across the torus, which we infer to 
exist close to the central engine. Conversely, the dust pos- 
tulated to be present in the NLR by their model is indeed 
detected with the AO data set. 

Granato & Danese (1994) and Granato et al. (1996, 
1997) developed a simple thick (te >30) torus model ex- 
tended over several hundreds pc. To minimize the number 
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of free model parameters, they have adopted a dust den- 
sity distribution constant with radial distance from the 
nuclear source. But, they do not rule out the possibility, 
in the case of smaller values of the optical depth (Te=1.5), 
that a more concentrated density distribution exist. Their 
predicted size and shape for the near infrared emission are 
compatible with those derived through AO observations at 
2.2, 3.5 and 4.8 /xm. 

A revised modeling of the AGN in NGC 1068 is timely, 
owing to the emergence of sub-arcsec resolution images in 
the near infrared (AO techniques) and in the millimeter 
range (interferometric techniques), giving direct access to 
the dust and molecular environment of the central engine. 

9. Conclusion 

The observation, for the first time at high angular reso- 
lution, of NGC 1068 at 3.5 and 4.8 /im provides new in- 
formations to build a more realistic model of this AGN 
under the current unification scheme. 
As regard to the AGN structure, we do observe: 
(i) an unresolved core, already known at 2.2 /xm to have 
a size (FWHM) less than 8 pc, and interpreted as the in- 
ner region of a dusty/molecular torus in which the central 
engine of NGC 1068 is embedded, 

(ii) along P.A.~100°, an extended emission up to 40 pc on 
either side of the core, particularly prominent at 4.8 ^m. 
This structure, also detected at 2.2 ^m up to 20 pc on ei- 
ther side of the core, is coincident in P. A. with the parsec- 
scale disc of ionized gas detected with VLBI (P.A.^110°), 
and is found to be roughly perpendicular to the axis of 
the ionization cone in NGC 1068. We interpret it as the 
trace, up to a 40 pc radius, of this dusty/molecular torus 
seen edge-on, 

(iii) an extended emission along the NS direction up to 
50 pc from the core and with rather symmetrical proper- 
ties on either side, both at 3.5 and 4.8 /xm. Again, this 
extended emission is detected , both at 2.2 /xm on a simi- 
lar scale, and at 10 and 20 /im on a slightly larger scale. 
It reveals the presence of dust in the NLR, heated both 
by hard radiation within the ionizing cone and by shocks 
associated with the AGN radio jet. 

As regard to the dust temperature and dust distribution 
in the central arcsec of the nucleus, we get a final pic- 
ture as follows. As close as r ~lpc from the central en- 
gine, the dust is heated up to its evaporation temperature, 
1500 K. Then the dust temperature declines very rapidly 
and reaches T'^500K at r < 28 pc. The total mass of 
warm dust within a 22 pc radius region is found to be 
around 0.5 M©. 

It is observed as well that the near infrared flux of 
NGC 1068, in the 2.2 to 4.8 /im range, has increased by a 
factor two over some 20 years, while in the 10 fim window 
the flux increase is only by a factor 1.2. 

Our results bring observational evidence of a dusty 
torus in the AGN of NGC 1068. They further support AGN 



modelling in the framework of the unification scheme: a 
thick torus surrounding a central engine (black hole and 
accretion disc). Although several models of the AGN in 
NGC 1068 are available, none matches in detail all the as- 
pects of the current near-infrared results obtained at a 
subarsec scale. Such new observational constraints make 
it both timely and exciting to run updated models. Yet, 
we are aware that the most convincing and undisputable 
argument for the presence of the torus-like structure will 
come from a study of the kinematics of the gas within the 
100 central parsec of NGC 1068. We expect such informa- 
tion to be soon obtained from ISAAC/ANTU observations 
on Paranal. 

Acknowledgements. We warmly thank J. P. Veran and E. Gen- 
dron for useful discussions, F. Lacombe for his help on the 
data reduction and Z. Tsetanov for his precious contribution 
in deriving the HST composite image seen by the WFAS. We 
acknowledge as well precious comments from an anonymous 
referee. 



References 

AUoin D., Marco O., 1997, Ap&SS, 248, 237 

AUoin D., Pantin E., Lagage P.O. et al., 1999, A&A, in prep. 

Antonucci R.R.J. , Miller J.S., 1985, ApJ, 297, 621 

Bailey J., Axon D.J., Hough J.H., Ward M.J., McLean I., et 

al., 1988, MNRAS, 234, 899 
Barvainis R., 1987, ApJ, 320, 537 

Beuzit J.L., Hubin N., Gendron E., et al., 1994, SPIE 2201 
Bland-Hawthorn J., Sokoloski J., Cecil G., 1991, ApJ, 375, 78 
Braatz J. A., Wilson A.S., Gezari et al., 1993, ApJ, 409, L5 
Bridger A., Wright G.S., Geballe T.R., 1994, in Infrared As- 
tronomy with Arrays: the next generation, McLean Eds., 

p. 537 
Capetti A., Macchetto F., Axon D.J., Sparks W.B., Boksen- 

berg A., 1995, ApJ, 452, L87 
Chapman S., Walker G., Morris S., 1999, in Proceedings of the 

ESO/OSA Sonthofen Topical Meeting on Astronomy with 

Adaptive Optics, 7-11 Sept. 1998, ed. D. Bonaccini, ESO, 

p. 73 
Chelli A., Perrier C, Cruz-Gonzales I., Carrasco L., 1987, 

A&A, 177, 51 
Clavel J., Wamsteker W., Glass I.S., 1989, ApJ, 337, 236 
Efstathiou A., Hough J.H., Young S., 1995, MNRAS, 277, 1134 
Efstathiou A., Rowan-Robinson M., 1995, MNRAS, 273, 649 
Evans I., Ford H., Kinney, A., et al., 1991, ApJ, 369, L27 
Gallimore J.F., Baum S.A., O'Dea C.P., Pedlar A., 1996a, ApJ, 

458, 136 
Gallimore J.F., Baum S.A., O'Dea C.P, Brinks E., Pedlar A., 

1996b, ApJ, 462, 740 
Glass I., 1991, MNRAS, 276, L65 
Glass I., 1997, Ap&SS, 248, 191 
Granato G.L., Danese L., 1994, MNRAS, 268, 235 
Granato G.L., Danese L., Franceschini A., 1996, ApJ, Lll 
Granato G.L., Danese L., Franceschini A., 1997, ApJ, 147 
Kriss G.A., Davidsen A.F., Blair W.P., Ferguson H.C., Long 

K.S., 1992, ApJ, 394, L37 
Krolik J., 1999, in Active Galactic Nuclei, Princeton Series in 

Astrophysics, Princeton University Press 



O. Marco & D. AUoin: Dusty/Molecular Torus in NGC 1068 



Lacombe F., Marco O., Geoffray H., et al., 1997, PASP, 110, 

1087 
Lynds R., Faber S.M., Groth E.J., et al., 1991, ApJ, 369, L31 
Marco O., AUoin D., Beuzit J.L., 1997, A&A, 320, 399 
Marco O., AUoin D., 1998, A&A, 336, 823 
Miller J.S., Goodrich R., Mathews W., 1991, ApJ, 378, 47 
Pension M.V., Pension M.J., Selmes R.A., Becklin E.E., 

Neugebauer G., 1974, MNRAS, 169, 357 
Pier E.A., Antonucci R., Hurt T., Kriss G., Krolik J.H., 1994, 

ApJ, 428, 124 
Pier E.A., Krolik J.H., 1992a, ApJ, 399, L23 
Pier E.A., Krolik J.H., 1992b, ApJ, 401, 99 
Pier E.A., Krolik J.H., 1993, ApJ, 418, 673 
Pogge R., 1988, ApJ, 328, 519 
Rieke G., Low F., 1975, ApJ, 199, L13 

Rouan D., Rigaut F., AUoin D., et al., 1998, A&A, 339, 687 
Thatte N., Quirrenbach A., Genzel R., Maiolino R., Tecza M., 

1997, ApJ 490, 238 
VeiUeux S., Goodrich R.W., HiU G.J., 1997, ApJ, 477, 631 
Young S., Hough J.H., Axon D.J., Bailey J.A., Ward M.J., 

1995, MNRAS, 272, 513 



